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Abstract 

We show that, though they are rare, there are asymptotically flat space- 
times that possess null geodesic congruences that are both asymptotically 
shear- free and twist-free (surface forming). In particular, we display the 
class of space-times that possess this property and demonstrate how these 
congruences can be found. A special case within this class are the Robinson- 
Trautman space-times. In addition, we show that in each case the congruence 
is isolated in the sense that there are no other neighboring congruences with 
this dual property. 
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I. Introduction 



For an arbitrary asymptotically flat space-time, in general, it is easy to 
find families of null geodesic congruences that are either twist-free or are 
asymptotically shear- free. To find congruences with both properties: i.e., a 
null geodesic congruences that is simultaneously twist-free (i.e., surface form- 
ing) and asymptotically shear-free is difficult and often impossible. By taking 
an arbitrary one parameter family of slicing of null infinity, 3 + , the inward 
directed null normals generate a twist-free null geodesic congruences. On the 
other hand, the asymptotically shear-free null geodesic congruences, (which 
are obtained from any shearing twist-free Bondi congruence via solutions of 
the good cut equation 1 ), almost always contain twist. It has proven diffi- 
cult to find examples of asymptotically flat space-times where one can find 
congruences that are both asymptotically shear-free and twist-free. They 
certainly are not generic. Special examples of such space-times are given 
by the class of regular Robinson- Trautman metrics, where the congruence is 
given by the twist-free degenerate principle null vectors of the Weyl tensor. 
It turns out that such space-times and their associated congruences are of 
considerable interest for studies of certain issues in the asymptotic behavior 
of gravitational fields. 1-3 It is the purpose of this note to show how large 
classes of such space-times can be constructed along with their associated 
asymptotically shear-free and twist-free congruence. In a future paper we 
will investigate their physical properties. 

II. Construction 

The basic problem is how to construct on 3 + ,(with given Bondi coordi- 
nates (u, C, £, t), special families of Bondi asymptotic shear, a(u, £, £), such 
that one can find from them a one parameter family of real solutions to the 
so-called 'good cut' equation, 

5 2 Z = a(Z,C,C). (1) 

In general, the solutions, Z(z a XX), °f the good cut equation depend on 
four complex constants, z a , {which define H-space}. For each arbitrarily 
chosen complex curve on the H-space, i.e., for z a = £ a (r), with complex pa- 
rameter r, one can easily construct an asymptotically shear-free null geodesic 
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congruence. 1 This is done by taking the solution, 



u^X(r,C,0 = Z(C(r),C,0 (2) 

and then defining (parametrically) the complex stereographic angle field, 
L(u,(,() on 3 + , by 

L(u,CX) = 5 (T) X(r,C,C) _ (3) 

u = X(t,(,0^t = T(u,(,(). (4) 

(5( T ) denotes the 5 operator holding r constant.) The shear-free congruence, 
defined and determined by the asymptotic vector l* a , is then obtained from 
the Bondi congruence by the null rotation, near 3 + , of their respective tangent 
vectors, (l* a ,l% ) by 



l* a = l a B -^m a B + ^-m a B + 0(r- 2 ), (5) 
m * a = n f B -L n ^ B + 0{r- 2 ), 

r?f a = m a B - -n a B + 0(r- 2 ), 
n* a = n% 

where (1%, n B , m B , m B ) form the Bondi null tetrad system at 3 + . The twist, 
S, of l* a is given by 

ill = -{6L + LL - 6L-LL }. (6) 

where dot denotes u-derivative. 

Our desire is to be able to chose real world-lines z a = £ a (s), r = s, so 
that the solutions to Eq.(l) take the form 

u = X(sX,C) (7) 

with both s and X real. Such solutions, if they could be found, define a 
(real) slicing (or NU coordinates) of 3 + , by s = const., where the null inward 
normals form null surfaces, i.e. are twist-free and are asymptotically shear- 
free. This later result follows from the fact that Eq.(7) is a solution to 
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the 'good cut equation'. 1 The difficulty has been that Eq.(l) is a complex 
differential equation and almost all solutions are 'intrinsically' complex. It 
was not at all clear what is the class of Bondi shears that allow for these 
one-parameter families of real solutions. 

The first thing to note is that for a completely general Bondi shear there 
are two real functions of (u, (, () that are buried in the (complex) shear, 
a(u, (, (), and we want one real function, u = X(s, (, (), to solve the equation 
- or from another point of view, 4 to 'transform the shear away'. We look for 
a Bondi shear that depends on only one real function. 

The construction of appropriate Bondi shears, though not immediately 
obvious, turns out to be quite simple. We begin by choosing an arbitrary 
real function of three variables, 

u = H(s,CX). (8) 

that can be interpreted as yielding a one-parameter slicing of 3 + . The inward 
directed null normals, as we mentioned earlier, form a null surface, i.e., are 
twist free by construction. Our task is to construct asymptotically flat space- 
times possessing a special class of Bondi shears so that the asymptotic shear 
of these normal null rays vanishes. 

We assume that Eq.(8) can be inverted as 

s = S(u,(,(). (9) 

Next (by applying 5 2 to H, holding s constant) we define the spin-wt-2 
complex 'electric type' function 

a*(s,CX) = ^ s )H(s,CO- (10) 

Finally the Bondi shear is obtained by replacing the s in Eq.(10 ) by s = 
S(u,(,() obtaining 

a(u,C,0 = a*(S(u,C,0,C,0 (11) 
so that the good cut equation becomes 

5 2 X = a(X,C,C) = tf s) H(S(XXX)XX)- (12) 

We easily see, from the construction and Eq.(10), that a real solution is 
given by 

u = X(s,(,() = H(s,CX)- (13) 
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We thus see that for any choice of the function u = H(s,(, C)j ( an NU 
coordinate transformation), we can define or construct a Bondi shear on 
the space-time so that the null normals to the constant 's' cuts of 3 + form a 
congruence that is both asymptotically shear-free and twist-free. The / = 0, 1 
harmonic components of H, written as £ a (s)/ a (C 0? determine the associated 
real world line, £ a (s). From the construction, it is clear that this, Eq.(ll), is 
the only class of Bondi shears with this dual property. 

Remark 1 In this construction it might appear that by simply adding to the 
H in Eq.( 13) another term of the form £* a (<s)^a(C C) a new solution could 
be constructed. This is not true since changing H by that term changes the 
S(u, C, C) an d thus the Bondi shear itself. The associated space-time would 
also change. However, if we are given the modified special shear from Eq.(ll) 

a(u, C, C) = (T*(S(u, C, C) , C, C) + d 2 a((, () (14) 

(which defines an arbitrary spin. wt. 2 function, F((,() = ( 5 2q; (C)C)) the extra 
term can be considered as pure gauge since the good cut equation, b^X = 
o-*(X,(,()+tfa((,(), becomes 

d 2 (X-a((,0 = d 2 X* = a*(X,(,0 

so that the gauge choice is equivalent tothe choice of a Bondi supertranslation. 

Another question that immediately arises: Given the Bondi shear, Eq.(ll), 
and the real one-parameter family of solutions, Eq.( 13), are there other real 
solutions close to the first real solution? 

For the generic situation the answer is no. 

To show this we begin with the good cut equation 

5 2 Z = a(Z,CC) 

and assume that the shear is of the form Eq.(ll) and that we are considering 
the real one parameter solution and the nearby solutions. We know from the 
general theory 5 ' 6 that there is a complex four parameter family of solutions, 
u = Z(z a , C, C) and hence the real one parameter solution has an associated 
real world-line, z a = £ a (s), with 

X(s,C, = Z(C(s), CO- (15) 
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For a given real point on the real world-line, z a = £ a (s ), we consider 
real neighboring points, z a = £ a (so) + <5£ a - More specifically we consider two 
different neighboring points; one on the world-line, 

z a = e(s )+e'(s )ds (i6) 

the other off the world-line, 

z a = e(s )+ V a (s )e (17) 

with 

V a (s )^^ a '(s ). (18) 

Given the solution, u = X(s ,CC) — Z(£ a (so), (, C) = Z , neighboring 
solutions, 

z(e+6e,co = z +d a z 5e 

are determined by the solutions, d a Z , to the linearized equation 

b\d a Z ) = &(Z ,C,0-d a Z . (19) 

The fact that we already know one solution, namely the one coming from 
displacement along the known world-line, Eq.(16), i.e., 

V = d a Z -e'(s ) = d s X(s,(X) (20) 

allows us the find the general solution. Prime denotes s-derivative. 
Assuming the general (four-parameter) solution has the form 

d a Z = VL a , (21) 

Eq.(19) becomes 

tf(VL a ) = 5V • L a + 2W ■ 5L a + V<S 2 L a = aVL a (22) 

or 

2WSL a + VtfL a = (23) 
5(F 2 5L a ) = 
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with solution 

V 2 <5L a = 7 >, (24) 

m = ^(i-f,-^(i + C 2 ), 2C) (25) 

Y a — three constant 'vectors' (26) 

Assuming (from genericity) that V~ 2 contains a non-vanishing I = 
spherical harmonic component, say K, we can write 

V' 2 = K + 1SW (27) 

with W a spin.wt.(-l) quantity, (containing, in general, all harmonics, I > 1). 
Using Eq.(27), Eq.(24) can be rewritten as 

= f a (K + bW) mi 

= jKKrrii + mfiW) 

= Y a (KNi + b[Wmi]) 

= ^(Kk + Wrm)} 

k = ^(C + C, -i(C-C),-i + CC) 

or 

L a = 1 ° a + f a (Kl l + Wm l ) 

with 7° another independent constant 'vector'. 
The general solution, i.e., Eq.(21), becomes 

d a Z = V{ 1 ° a + fa(Kk + Wm l )}. (28) 

The four constant 'vectors', (7°, 7a), are determined (partially) by the con- 
dition that, 

V = d a Z,C' = V{C'll + C'Y a (Kh + Wm t )} (29) 

or 

C'la + C'YaiKk + Wm,)}^!. 
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Since (l^W^mj) have angular behavior the four 'vectors' must satisfy the 
algebraic conditions 

C'7° = 1 (30) 

C'Ya = 0. (31) 

Assuming that we have found a set of four independent 'vectors' with 
the given £ a/ that satisfy Eqs.(30) and (31), then for any other real quantity, 
rj a (sq) , with 

V a (so) ^ rf a \s ) (32) 

we have that 

V a (s )f a ? 0. 

This immediately implies that all solutions that differ from V or a multiple 
of V, i.e., 

d a Z oV a = y{rf 7 ° + V a fa(Kh + Wm t )} 

are intrinsically complex, so that in the generic situation there are no other 
real solutions in the neighborhood of Eq.(13). 

III. Conclusion 

We studied here the question of: what is the class of space-times that 
possess null geodesic congruences that are both asymptotically shear-free 
and twist-free and furthermore asked, within that class how common are 
such congruences? 

Our motivation for this work was the recognition that the well-known 
Robinson-Trautman metrics allow such a congruence - and in fact are essen- 
tially defined by a special case in this class, namely the ones that are globally 
shear-free and twist-free and not just asymptotically shear-free and twist-free. 
This work can thus be considered as attempting to find a generalization of 
the Robinson-Trautman case. 

In a future paper we will apply the results of this work to study the 
relationship of the real world-line to the Bondi energy-momentum four- 
vector and the resulting 'equations of motion' 3 coming from the asymp- 
totic Einstein-Maxwell equations, i.e., those coming from the Bondi energy- 
momentum loss equation. Though it might be hard to believe, we are con- 
vinced that this future study can and will shed light on the run-away behavior 
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associated with the radiation reaction term in the electrodynamic equations 
of motion for a charged particle. 3 
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